Chaotic dynamics of a Bose-Einstein condensate coupled to a qubit 
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We study numerically the coupling between a qubit and a Bose-Einstein condensate (BEC) moving 
in a kicked optical lattice, using Gross-Pitaevskii equation. In the regime where the BEC size is 
smaller than the lattice period, the chaotic dynamics of the BEC is effectively controlled by the qubit 
state. The feedback effects of the nonlinear chaotic BEC dynamics preserve the coherence and purity 
of the qubit in the regime of strong BEC nonlinearity. This gives an example of an exponentially 
sensitive control over a macroscopic state by internal qubit states. At weak nonlinearity quantum 
chaos leads to rapid dynamical decoherence of the qubit. The realization of such coupled systems 
is within reach of current experimental techniques. 
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The dynamics of Bose-Einstein condensates (BEC) in 
the regime of quantum chaos is now actively investigated 
by several experimental groups in the world [l], 0, S 0] • 
These experiments implement the quantum version of 
the Chirikov standard map (kicked rotator) [B| by means 
of kicked optical lattices. This model exhibits interest- 
ing phenomena such as dynamical localization of diffu- 
sive energy growth, which has been first observed in cold 
atoms [6]. However, for cold atoms the dynamics is de- 
scribed by the linear Schrodinger equation, while for BEC 
nonlinear effects are present that can be studied in a 
wide range of situations with Gross-Pitaevskii equation 
(GPE) 7] . Due to this nonlinearity the dynamics of BEC 
can become truly chaotic with appearance of exponential 
instability [1, Q , provided the size of BEC is smaller than 
the optical lattice period. In many cases, a correct de- 
scription of the system is obtained by accounting for the 
internal structure of the atoms forming the BEC. Such 
BEC must be seen as quantum objects with several com- 
ponents. In the situation where there are two compo- 
nents, the system becomes equivalent to a BEC coupled 
to a qubit. This leads to the fundamental problem of the 
dynamics of a truly chaotic BEC coupled to a qubit and 
controlled by it. At the same time, this dynamics gener- 
ates a feedback on the purity and evolution of the qubit, 
which may significantly affect its properties. Since the 
BEC dynamics is exponentially unstable, it has the po- 
tential to become a detector with exponential sensitivity 
to the qubit state. The spinor states of BEC are now ac- 
tively studied experimentally (see e.g. 1(J U|), such that 
in principle it would allow to realize the coupling between 
a qubit and chaotic BEC dynamics. The decoherence of 
a qubit coupled to a BEC has also been discussed in the 
context of solid state physics devices [l2j |. 



To describe the dynamics of a BEC coupled to a qubit 



in a kicked optical lattice we use GPE Q 
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The first term on the r.h.s. accounts for the free evo- 
lution of the cold atoms, the second one represents the 
effect of the optical lattice of period A = 2it/k, 5t(X) 
being a periodic delta-function with period T, and the 
third one corresponds to the driving of the qubit with a 
strength S. The last term describes the nonlinear inter- 
actions between atoms treated in the framework of GPE 
for spinor BEC. We start considerations with the case 
gu = gio — goi = goo = g, where the nonlinear parame- 
ter g — NgiD is determined by the number N of atoms 
in the condensate and the effective ID coupling constant 
giD — —2aohiu±, with lo±_ the radial trap frequency and 
ao the 3D scattering length. The spinor wavefunction is 
normalized to one (L(|'0o( a O| 2 + \^pi{x)\ 2 )dx = 1). In the 
following, we express the momentum of atoms in recoil 
units and time t in units of the kick period T by setting 
H = m = k = 1 in Eq. |T|). In absence of atom-atom 
interactions at g = 0, Eq. (Q]) reduces to the usual kicked 
rotator model with classical chaos parameter K = kT 
and effective Planck constant T (see e.g. Q). The 
first three terms on the r.h.s. correspond to the linear 
evolution of the kicked rotator coupled to a qubit. The 
strength of the kick depends on the orientation of the 
qubit due for example to different detunings between the 
atomic transition frequencies and the laser frequency 0] . 
This linear model was studied in [pj ] where it was shown 
that quantum chaos in the kicked rotator leads to rapid 
decoherence of the qubit. Here we study the effects of a 
nonlinearity g induced by interactions between atoms in 
the BEC. 

We consider an initial state of the form 1^) = \ib) ® \<h) 
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FIG. 1: (Color online) Probability density \^i(x)\ 2 at t = 
10 for k = 3.4, e = 0.5, T = 2, 5 = 0.2 showing a broad 
distribution at g — (red/grey curve, magnified by a factor 
100) and a BEC soliton at g = 20 (blue/black curve). The 
initial state is given by Eq. ((2)1 at t — for go = 20, xo = 0.1 
andpo = 0.05 and with qubit state \<j>) = (|0) + \l))/y/2. Inset 
shows data on a smaller scale near the soliton center, with full 
curves for |?/ii(a;)| 2 and dashed curves for |f/>o(£)| 2 . 



where = a\0) + j3\l) is the internal state of the BEC 
and its motional state given by the soliton distribu- 
tion 

w ^ = y/go exp (ipojx -xq- ppt/2) + zggt/8) 
2 cosh(g (x - x -p t)/2) 

taken at time t — 0. For k — and g = 9 this is the 
exact solution of Eq. |T]), which describes the propagation 
of a soliton with constant velocity po 



14| . Numerical 



simulations are performed using the second order SABA 
integrator [15J for which the error scales as 0(At 3 ) with 
time step At. Computations have been done with Ap = 
1/64, total number of states N$ — 2 15 — 2 16 and time 
step At « 0.01/(3+1). 

At moderate values of k, a narrow distribution with a 
solitonic form ([2]) spreads rapidly in space and momen- 
tum in absence of interactions between atoms (g = 0). 
In contrast, for sufficiently large g, the shape of the soli- 
ton is only slightly perturbed even after a large number 
of kicks (see Fig. [T|). In this regime, the center (x c ,p c ) 
of the BEC follows the dynamics described by Chirikov 
standard map p c = p c + ksmx c ; x c — x c +p c T, 
where bars denote the values of the soliton position and 
velocity after one kick iteration. Accordingly, the soliton 
dynamics becomes chaotic for K = kT > 1. In presence 
of the qubit dynamics imposed by the Rabi splitting 5, 
the kick function f{xt) — (%t+i ~ % x t + xt-i) / (kT) starts 
to depend on the internal qubit state (here Xt is the av- 
erage position of all atoms from the two components at 
time t) . This leads to a splitting of / amplitude between 
the two values (1 ± e/fc)sina; determined by the qubit 



FIG. 2: (Color online) Circles show the kick function f(x) 
obtained from the BEC dynamics of Eq. (TfJ) with parameters 
k = 1, e = 0.25, T = 4, 8 = 0.2 and g = 20. Initial state is 
the same as in Fig. [T] The black solid line shows f(x) = sin a; 
while the grey area, induced by the qubit, corresponds to the 
region delimited by the curves (l±e/fc) sin x. Inset: Poincare 
section of the standard map at K = 1 (small dots), BEC 
positions are shown by large red/grey dots for parameters 
k = 1, e = 0.5, T = 1, 5 = 0.05 and g = 20 and initial state 
as in the main Figure. 



state. In presence of Rabi splitting <5, the amplitude / 
varies between these two limiting values due to nontrivial 
qubit dynamics as is clearly seen in Fig. [51 Due to the 
qubit oscillations, the soliton can even penetrate inside 
classically integrable regions (see inset of Fig. [2]). Inside 
the chaotic region a soliton trajectory in the phase space 
is exponentially sensitive to the qubit internal quantum 
state. 

The GPE nonlinearity g also influences strongly the 
qubit dynamics. Indeed, for g = the purity of the 
qubit £(t) = |Tr (p(t)a) | decays rapidly to almost zero 
due to quantum chaos in the dynamics of atoms, as can 
be seen in Fig. [3] The decay is linked to the Lyapunov 
exponent of classical chaotic dynamics as discussed in 
[1 31 ] - In contrast, for large values of 3, the purity of the 
qubit decays very slowly, remaining close to 1 after al- 
most one hundred kicks (see Fig. [3]). In this regime, the 
purity decay is determined by the decay of the solitonic 
part of the BEC. Indeed, the latter can be defined as 
W ' = C-w/g(\Mx)\ 2 + \M*)\ 2 )dx with x c taken as 
the soliton center position, and the data in Fig. [3] show 
that it decays in the same way as the purity £. Since the 
lifetime of the soliton t s can be rather large (predicted 
to be t s ~ g 4 /k 2 in 8]), this implies that the qubit can 
remain pure for a very long time. 

In the regime where the qubit purity is high, the 
probability oscillates between the two internal states of 
the BEC. This leads to oscillations of the polarization 
of the qubit £ z = Tr (pcr z ), which are clearly seen in 




FIG. 3: (Color online) Qubit polarization £ z (top and middle 
panels) and purity £ (bottom panel) as a function of time for 
k = 3.4, e = 0.5, T — 2, and 8 = 0.2. Data are shown for: (top 
panel) gn = gio = goi = goo = g with g — 20 (red triangles) 
and 3 = (blue circles); (middle panel) goo = gn = 21 and 
9oi = gio = 19 (orange diamonds) and goo = gn = 30 and 
goi = <?io = 10 (green squares). In the bottom panel, the 
black solid line shows the solitonic part W s (see text), other 
symbols are as in top and middle panels. Initial state is the 
same as in Fig. [1] 
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FIG. 4: (Color online) Distance 8x between solitons with the 
qubit in the states up ipi(x,t) and down tpo(x,t) as a func- 
tion of time t for goo = gn = g = 40, goi = gio = (red 
squares), goo = ffn = g = 39.5, goi = gio = 0.5 (green cir- 
cles), goo = gn = g = 28, goi = gio = 12 (blue diamonds), 
and goo = gn = goi = gio = g = 20 (orange triangles). Other 
parameters are as in main part of Fig. [T] 



Fig. [3] These oscillations remain stable with respect 
to small variations of coupling matrix elements gij (e.g. 
.9io = .9oi = 19 < ,9oo = .9ii = .9 = 21 in Fig. [3]). How- 
ever for a significant asymmetry 510 = .901 = 10, 300 = 
.9ii — .9 — 30 the nonlinear shift strongly modifies the 
unperturbed Rabi frequency and the amplitude of oscil- 
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FIG. 5: (Color online) Density plot of the spectral density 
S(u) of £ z (t) expressed in arbitrary units and shown as a 
function of the frequency v and the nonlinear parameter g 
for k = 3.4, s = 0.5, T = 2, 5 = 0.2, and the rescaled Rabi 
frequency vr = T5/tc. Initial state is the same as in Fig. [T] 

lations of £, z (t) becomes strongly suppressed, even if the 
purity remains close to unity. In this regime the distance 
between the centers of the two solitons with the qubit 
in up and down states remains small compared to their 
size (see Fig. [4|. Only for small values of gtj = g the 
purity starts to drop rapidly and the oscillations of £ z (t) 
decay quickly due to dynamical decoherence induced by 
quantum chaos. When the coupling matrix has small 
off-diagonal elements (e.g. g i = gio < 0.5, see Fig. H]), 
the solitons for the two components move independently 
from each other and the purity drops to zero as soon as 
the solitons are separated. 

To better analyze the properties of these qubit oscil- 
lations, we determined the variation of spectral density 
S(v) of £ z (i) as a function of g. This spectrum is shown 
in Fig. [5] For small values of g, the spectral density is 
very low since quantum chaos leads to a rapid decay of 
qubit oscillations and purity. In contrast, for larger val- 
ues of g > 10 a prominent peak at the Rabi frequency 
Vji = T5/tt becomes well pronounced, even if the soliton 
dynamics is strongly chaotic (K = kT = 6.8). 

The interplay of quantum chaos and GPE nonlinearity 
and their effects on the purity of the qubit are displayed 
in Fig. [5] For small g < 10, the purity is destroyed 
with the increase of the chaos strength represented by 
the parameter K. On the contrary, for g > 10 the purity 
remains close to 1 even for large values of K when the 
soliton dynamics is strongly chaotic. 

This behaviour can be qualitatively understood by 
the following arguments. From Eqs. (P)-© the cou- 
pling energy of the single-component soliton at go = g 
is E s = -g 2 /& [16]. For g n « g w « g Q1 w g 00 w g > 1 
the energy of attraction between the two components 
is of the same order of magnitude, thus being signifi- 
cantly larger than the difference E). between the kinetic 
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FIG. 6: (Color online) Density plot of the qubit purity £ after 
t = 20 kicks as a function of the chaos parameter K = kT 
and the nonlinear parameter g for T = 2, S — 0.2 and e = 0.5. 
Initial state is the same as in Fig. [T] 



energies induced by the kick term in Eq. ([T]). Indeed, 
the latter can be estimated as Ek ~ pAp ~ 2ke\/i; for 
fc « 3, t = 20, e = 0,5 this gives Ek ~ 13 and the con- 
dition E s « Ek gives the critical value of g c « 10. This 
estimate for g c is compatible with the numerical results 
displayed in Fig. [6l even if this approach is rather sim- 
plified. In addition, we should also note that according 
to the well-known result 14J the initial state ^ has a 
coupled solitonic state for g > 0.18<?o- F° r the conditions 
of Fig. [6] this gives a critical g w 3.6 that is significantly 
smaller than the numerical value of g c ~ 10. We at- 
tribute this difference to the internal dynamics of qubit 
which is not taken into account by the estimate above. 
In the case g\\ = goo = g ^> goi = <?io the coupling 
energy between solitons can be obtained by degenerate 
perturbation theory which gives E s ~ — <7oi<7- I n the case 
where Ek > E s , the solitons start to separate exponen- 
tially fast due to classical chaotic dynamics. Such kind 
of behaviour can be potentially used to develop detectors 
with exponential sensitivity to the qubit state. 

In conclusion, our results show that the coupling of 
a qubit with a chaotic BEC dynamics leads to rather 
nontrivial behaviour. On one hand, the qubit can con- 
trol the dynamics of the BEC, moving it from chaotic 
to integrable regimes and vice versa. The BEC position 
is exponentially sensitive to the state of the qubit. In 
spite of chaotic motion of BEC, the purity of the qubit 
is well-preserved in this situation. Such an exponential 
sensitivity of a global BEC position to an internal quan- 
tum state appears due to the absence of second quan- 
tization for BEC motion in GPE. In presence of such 
quantization the whole system is described by a linear 
Schrodinger equation and the chaotic dynamics of BEC 



can continue only on a relatively short Ehrenfest time 
scale as it is discussed in [9]. On the other hand, in the 
regime of weak BEC nonlinearity, quantum chaos leads 
to rapid dynamical decoherence of the qubit. The real- 
ization of such coupling between a qubit and a BEC with 
chaotic dynamics is within reach of current experimental 
techniques. 
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